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The ISW-tSZ cross correlation: ISW extraction out of pure CMB 
data 
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ABSTRACT 

If Dark Energy introduces an acceleration in the universal expansion then large scale 
gravitational potential wells should be shrinking, causing a blueshift in the CMB photons that 
cross such structures (Integrated Sachs-Wolfe effect, [ISW]). Galaxy clusters are known to 
probe those potential wells. In these objects, CMB photons also experience inverse Compton 
scattering off the hot electrons of the intra-cluster medium, and this results in a distortion 
with a characteristic spectral signature of the CMB spectrum (the so-called thermal Sunyaev- 
Zel'dovich effect, [tSZ]). Since both the ISW and the tSZ effects take place in the same poten- 
tial wells, they must be spatially correlated. We present how this cross ISW-tSZ signal can be 
detected in a CMB-data contained way by using the frequency dependence of the tSZ effect 
in multi frequency CMB experiments like Planck, without requiring the use of external large 
scale structure tracers data. We find that by masking low redshift clusters, the shot noise level 
decreases significantly, boosting the signal to noise ratio of the ISW-tSZ cross correlation. 
We also find that galactic and extragalactic dust residuals must be kept at or below the level of 
~ 0.04 (yuK) 2 at I — 10, a limit that is a factor of a few below Planck's expectations for fore- 
ground subtraction. If this is achieved, CMB observations of the ISW-tSZ cross correlation 
should also provide an independent probe for the existence of Dark Energy and the amplitude 
of density perturbations. 
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1 INTRODUCTION 

The primary Cosmic Microwave Background (CMB) and espe- 
cially its angular power spectrum provides us with powerful con- 
straints on the content of the universe and its evolution. It is now 
well established that an accurate understanding of the primary 
CMB power spectrum requires a good understanding of the sec- 
ondary CMB anisotropies resulting from the interaction of the 
CMB photons with the matter along the l ine of sight from the 
last scattering surface to the observer (see I Aghanim et all I2008L 
for a review). The great efforts undertaken to understand these sec- 
ondary anisotropies, in order to best recover the primary CMB, also 
provide us with powerful independent cosmological probes when 
the secondary anisotropies are regarded as a source of information 
rather than contamination. 

Among those secondary CMB anisotropies, some result from 
the gravitational interaction of the CMB photons with the po- 
tential wells they cross. One of them is the Integrated Sachs- 
Wolfe (ISW) effect, by which CMB photons experience some 



E-mail: nicolas.taburet@ias.u-psud.fr 



blue/redsh ift as they pass throug h large scale time evolving poten- 



tial wells dSachs & Wolfd[l967 ). Since a dark energy like compo- 
nent is expected to affect the growth of large scale structures, mak- 
ing them shallower, a detection of the ISW effect is an important 
probe for establishing its existence - provided that the Universe is 
flat and general relativity is a correct description of gravity - and 
constraining the equation of state of such a component. 

Detection claims of the ISW effect arose as soon as the 
first year WMAP data were released. Those claims were based 
upon cross correlation analyses of WMAP CMB data and galaxy 
density templates built from different surveys. While most of 
the first analyses were conducted in real space (i.e., by comput- 
ing the angular cross-correlation function), subsequently new re- 
sults based upon Fourier/multipole and wavelet space were pre- 
sented. The results from the WMAP team on the cross corre- 
lation ofJ^R^VeryLarge Sky Survey (NVSS) with WMAP 
data dNolta et alj 120041) were soon followed by other analy- 
ses applied not only on NVSS data, but also on X-ray and 
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sequent data releases from both the CMB and the SDSS side be- 
came public, new studies prompted further evidence for signifi- 
cant cross correl ation between CMB and Large Scale S t ructur e 
(LSS) data, e.g.. IPadmanabhan et all j2005h: ICabre et al.l d2006h : 
iGiarmantonio et al.l d2006h : iRassat et all d2007l) . By that time, 
wavelet techniques were also applied on NVSS and WMAP data, 
providing the highest significance ISW detection claims at the level 
of 3- 4 cr, dVielva et all2006l:|Pietrobon et al.l | 2006l : lMcEwen et al.1 
|2007|) . The initial effort of lGranett et al.l d2008n . consisting in stack- 
ing voids and superclusters extracted from SDSS data, yielded a 
very high signifi cance (~ 4cr) ISW detection claim. However, it 
was later found in lGranettetailj2009h that such signal could not be 
due to ISW only, since a gravitational potential reconstruction from 
the Luminous Red Gala xy (LRG) sample of S PSS y ielde d a much 
lower signal (~ 2.5cr). iGiannantonio et al] d2008l) and iHo et al.l 
d2008h used different LSS surveys in a combined cross-correlation 
analysis with CMB data and claimed high significance (~ 4 — 5<r) 
ISW detections. 

However, doubts on the validity of such claim s have also 
arised recently. iHernandez-M onteagud o et al] d2006al) first pointed 
out the lack of significant cross correlation between WMAP 1st 
year data and density surveys built upon 2MASS, SDSS and NVSS 
surveys on the large angular scales, but detected the presence 
of radio point source emission and thermal Sunyaev-Zerdovi ch 
effect on the small scales. In Hernandez-Monteagudo (2008) a 
study of the expected signal-to-noise ratio (S/N) for different sky 
coverages was presented, and it was found that in the standard 
ACDM scenario the ISW - density cross correlation should be 
well contained in the largest angular scales, (I < 50 - 60). This 
was proposed as a consistency check for ISW detection agains t 
point-source contamination. In iHernandez-Monteagudol (2009) 
cross-correlation analyses between NVSS and WMAP 5th year 
data provided no evidence for cross-correlation in the large angular 
range (/ < 60). A signal at the 2-3 a level was however found at 
smaller scales, although its significance increased with increasing 
flux thresholds applied on NVSS sources (in contradiction with 
expectations for the ISW probed by NVSS and raising the issue of 
radio point source contamination). Furthermore, the intrinsic clus- 
tering of NVSS sources on the large scales (relevant for the ISW) 
was found too high for the commonly assumed reds hift distribution 
for NVSS sources, as found in previous works, dNegrello et al] 
120061 : iRaccanelli et al.l l2008h . Regarding ISW detection claims 
based upon SDSS data, there is also some ongoing discussion after 
recent failures in finding any statistical significance for the ISW, 
dBielbv et alj20icl : lLopez-Corre doira et al.l201ol ; ISawangwit et al] 
l20ld) . This situation is partially caused by the fact that the ISW 
is generated on the large angular scales and at moderate to high 
redshifts (z e [0.1,1.3]). Deep galaxy surveys covering large 
fractions of the sky are hence required to sample the ISW properly, 
but those are not available yet (or not properly understood). 

Ideally one would try to find the ISW contribution to CMB 
anisotropies by using CMB data exclusively. And in this context the 
thermal Sunyaev-Zel'dovich effect be comes of relevance. The ther - 
mal Sunyaev-Zel'dovich (tSZ) effect dSunvaev &Zel'dovichll972h 
results from the inverse Compton scattering of the CMB photons 
off the galaxy clusters electrons, and is expected to provide cosmo- 
logical constraints on the normalisation at SIT 1 Mpc of the den- 
sity fluctuations power spectrum, erg, as well as on the amount 
of matter Q.„, an d to a lower extent o n the dark energy equa- 
tion of state (e.g. iBattve & Welledl2003l) . Since both the tSZ and 
the ISW effects probe large scale structures and their evolution, 



a correlation is therefore expected between these two signals. As 
Hernandez-Monteagudo & Sunvaevl d2005l) pointed out, provided 
that the tSZ has a definite and well known frequency dependence, 
it is possible to combine different CMB maps obtained at different 
frequencies in search for a. frequency dependent ISW-tSZ cross- 
correlation. The advantages of this approach are two folded: (i) 
only CMB data (obtained at different frequency channels) are re- 
quired, and hence there is no need for using and characterising an 
external LSS catalog, and (ii) a better handle on systematics is pro- 
vided since the ISW-tSZ cross correlation has a perfectly known 
frequency dependence that can be searched for in multiple channel 
combinations. This approach in experiments like Wane^lJ, covering 
the whole sky in a wide frequency range, is well suited to separate 
the tSZ from other components present in the microwave range. 

After introducing in Section 2 the theoretical background, we 
study in Section 3 a combination of CMB maps at different fre- 
quencies that provides an unbiased estimate of the ISW-tSZ angular 
power spectrum. We also compute the expected significance for the 
ISW-tSZ cross correlation. In Section 4 we analyze how the large 
scale structure contributes to the tSZ and ISW auto power spectra 
in different redshift ranges. This allows us to define a strategy to 
optimize the S/N of the ISW-tSZ cross-correlation by applying a 
selective mask on galaxy clusters. We then investigate in Section 5 
the limitations of our method due to the presence of galactic and ex- 
tragalactic foregrounds and suggest some approaches to minimize 
their impact. We present our conclusions in Section 6. 



2 ANGULAR POWER SPECTRA OF LARGE SCALE 
STRUCTURES TRACERS 

CMB photons can interact with the matter situated along the line of 
sight from the last scattering surface to the observer, and thus suffer 
gravitational or scattering effects. This produces fluctuations of the 
observed CMB temperature on the sky. The most important scat- 
tering effect is known as the Sunyaev-Zel'dovich effect that arises 
when CMB photons scattered off the electrons of the intra cluster 
gas. Since these electrons are lying within the large scale structure 
potential wells, a correlation is thus expected between the SZ effect 
and the ISW temperature fluctuations due to the energy change of 
the CMB photons that pass through time evolving potential wells. 



2.1 Late ISW 

In a ACDM scenario, the accelerated expansion of the universe 
makes large scale potential wells to shrink. As a consequence, 
CMB photons crossing a potential well cf> do not loose as much 
energy exiting this well as what they gained when falling into it. 
This is known as the late ISW effect and results in a modification 
of the CMB blackbody temperature (e.g.. lMartinez-Gonzalez et"aT1 



AT 



2 r d^ 

--J d "dv 



(1) 



where 77 is the conformal time defined as dn = dt/a(t), where t 
and a{i) are the coordinate time and the scale factor in a standard 
Friedmann-Robertson- Walker (FRW) metric. The relation between 
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the gravitational potential and matter distribution variations is given 
by the Poisson equation in physical coordinates : 



V 2 </> = 4nGp m (l +6). 



(2) 



This equation is easier to solve in the comoving frame and in 
Fourier space (k stands for the wavenumber and the subscript c 
denotes comoving units): 
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(3) 



where H is the Hubble constant. We have used the relation 
between the critical density of the Universe and the expansion 
rate of the homogeneous background p c (a) = (3H 2 (a))/(8nG) = 

p m (a)l Ulmocr 3 (gjgj) ). 

In the linear regime, as long as the different k modes are not 
coupled to each other, the matter overdensity for a pressureless fluid 
can be written 8k{a) = D + (a)St(z = 0) where the growing mode 
D + (a) is a solution to the following differential equation : 



8(k) + 2HS(k) = 4nGp m 8(k), 



(4) 



in which dots denote a derivation with respect to the physical time. 

Using Eq. (0 we can express the CMB temperature fluctua- 
tions due to the ISW in the linear regime : 
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The temperature fluctuations due to the ISW effect can be projected 
on the spherical harmonics basis, 
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with the ISW transfer function 

J2 



(6) 
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2.2 The thermal SZ effect 

The SZ effect consists of two terms. The main one is the thermal 
SZ effect (tSZ) that is due to the inverse Compton scattering of the 
CMB photons off the intracluster gas hot electrons. The second one 
is the kinetic SZ effect (kSZ) which is a Doppler shift due to galaxy 
clusters' motion with respect to the CMB rest frame. The tSZ effect 
transfers some energy from the hot electrons to the CMB photons. 
As a result, in the direction of the cluster, the CMB intensity is de- 
creased in the Rayleigh Jeans part of the spectrum and increased in 
the Wien part. This translates into a characterisitic spectral signa- 
ture g v of the induced CMB temperature secondary fluctuations, ex- 
pressed as a function of the adimensional frequency x = hv/(kBT c ) 
and the comptonisation parameter y. Neglecting relativistic correc- 
tions, this parameter can be written as 
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The comptonization parameter y corresponds to the integrated elec- 
tronic pressure along a given line of sight through the cluster : 



dln.T, 



(9) 



The kSZ, on the other hand, does not have a different spectral signa- 
ture from the CMB. The secondary temperature fluctuations power 
spectrum due to the kSZ is about 2 to 4 times smaller than the one 



induc ed by the tSZ at 150 GHz l lSehgal et alj|2010l : iLueker et al.1 
2010). We can safely neglect its contribution in our analysis as ex- 
plained in Section 3. 

In order to calculate the temperature fluctuations due to 
a population of N c l usters, one can use t he halo approach 
dCole & Raised 1 19881 ; ICoorav & Shetr] l2002h . In this article, 
we adopt the line of sight approa ch that was introduced in 
iHernan dez-Monteagud o et alj |2006b). 
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where r p stands for physical distances, w for the line of sight posi- 
tion, and ym, is the value of ^rdTii, at the center of the cluster 
and u is the electronic radial pressure profile, >3d(>v) = y3D,c u ( w )- 
Replacing the discrete summation with an integral over the position 
we get : 
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The convolution between the mass function and the cluster profile 
is easier to handle in Fourier space : 
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Since galaxy clusters are not exclusively Poisson distributed, in the 
Fourier description of the spatial distribution of these sources a cor- 
relation term adds to the Poisson term. It represents the modulation 
of the cluster number density by the underlying density field : 
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dM 



dn 

(k D , M, r„) = 

, p p dM 
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where the l inear bias is mode lled with b(M,z) 
S c /cr 2 (R,z) dMo&WhiteMl996l) . 

That is, the power spectrum of the tSZ angular anisotropy 
can be written as the cont ribution of two different terms, (e.g., 
Komatsu & Kitavamj [l999). The first one refers to the Pois- 
son/discrete nature of these sources, and is known as the 1-halo 
term. The second one refers to the spatial modulation of the den- 
sity of these sources, obeying the large scale density field, and is 
referred to as the 2-halo term, which is sensitive to the underlying 
matter density distribution. In this article we are mainly interested 
in the detection of the cross ISW-tSZ term, to which only this sec- 
ond term contributes. 

The temperature fluctuations due to the SZ effect 2halo term 
can be projected on the spherical harmonics basis, 
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with the tSZ 2-halo term transfer function : 

dn 
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D + (r p )jf(k p r p ) 
The Fourier transform of y3o(r) is : 

y 3D (k p , M, z) = 4n J y 3D (r p , M, z)jo(k p r p y-dr p (16) 

Switching to comoving units, we write the a (m in a form similar to 
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Eqs. d6} and Q 



■° J 

■m transfer : 
g, f air. f ± 



r (m (KK Z2h (k c )s kc , z=0 



(17) 



with the tSZ 2-halo term transfer function : 



(M, r c )fc(Af, r c ) 



y 3D (fc P , M, r c )Z> + (>"c)./XVc) (18) 

In the f ollowing we have used th e ISheth et alj J200lh mass func- 
tion, the lKomatsu & Seliakl d2002l) model for the intracluster elec- 
tronic distribution and the WMAP-5 year cosmological parameters 
dKomatsu et alJl2009l) : fi A = 0.721, n m = 0.279, fi b = 0.046, 
h = 0.701, o-g = 0.817 and n s = 0.96. 




0.010 
0.001 



2.3 ISW-tSZ Cross power spectrum 

Using Eqs. JSJ, Q, J17t and J 1 3b it is straightforward to obtain the 
cross correlation between the ISW effect and the tSZ effect : 



2 J ZMeAfVfh 



)Af Z2h (k c )P(k c ,z = 0) 



where we have introduced the matter power spectrum : 
(s K:=0 S' k , J) = (2n) 3 S D (k-k')P(k,z = 0) 



(19) 
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We speed up the computation of the power spectra by adopting the 
Limber approximation : 
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This Limber approximation is based on the closing relation of the 
spherical Bessel function, 



k dkj e (kri)j e (kr 2 )T(k) =s = T\ , (22) 

where T is a function that slowly varies with k. We used k = ^-^2 
in order to ensure an error in 0(i~ 2 ) while the frequently em- 
ployed flat sky approximati on k = - only ensures an error in 0{l~ l ) 
i lLoverde & Afshordill2008l) . 

In Fig. [j]we present the CMB-CMB, ISW-ISW, tSZ-tSZ and 
ISW-tSZ power spectra calculated at 100 GHz without using the 
Limber approximation. We also represent the instrumental noise for 
a Planck like experiment as well as the cosmic variance associated 
totheCMB. 

Fig. [2] represents the relative difference between the various 
angular power spectra calculated using the Limber or the exact for- 
mula. While the 'standard' Limber approximation (k = Ijr) en- 
sures a 10% precision from t > 6 and i > 20 for the tSZ-tSZ 
and ISW-ISW power spectra respectively, these limit multipoles 
become I > 3 and t > 1 when using the k ^ (£ + l/2)/r approx- 
imation. As discussed in Afshordi et al] d2004l) the errors induced 
by this later approximation are small compared to the cosmic vari- 
ance. Therefore in the following we use this k = (£ + 1 12) I r version 
of the Limber approximation. 

The relative amplitude of the ISW-tSZ power spectrum 
with respect to the contaminant spectra (primary CMB and tSZ) 



Figure 1. ISW-tSZ angular cross power spectrum at 100 GHz (thick solid 
red line), ISW angular power spectrum (thin solid blue), tSZ-lhalo term 
(short dashed green), tSZ-2halos term (long dashed green). The CMB power 
spectrum is represented as a black dotted line for comparison. The dot- 
dashed line represents the cosmic variance associated to the CMB and the 
triple dot-dashed line is the Planck noise at 100 GHz. 




Figure 2. Relative error when using the k a (t + l/2)/r Limber approxi- 
mation (thick lines) or the k a tjr Limber approximation (thin lines). The 
solid red line stands for the ISW-tSZ angular cross power spectrum, the 
triple dot-dashed blue line for the ISW angular power spectrum, the dot- 
dashed green line for the tSZ-lhalo term and the dotted green line for the 
tSZ-2halos term. This plot clearly illustrates the great advantage of using 
the k ^ (t + l/2)/r approximation instead of the 'standard' k a t jr Limber 
approximation in low t studies of the power spectra. 



makes it less difficult to be detected at large scales t han a t 
smaller scales (more details in Section [3] s ee also lCooravl d2002l) ; 
iHernandez-Monteagudo & Sunva ev (2005)). For multipoles lower 
than 100, we see in Fig. [Tj that the effect of the instrumental noise 
can be neglected since its power spectrum is more than one order 
of magnitude lower than the astrophysical signals considered here. 
It is nevertheless essential to notice that the ISW-tSZ signal, we are 
interested in measuring, is much smaller than the cosmic variance 
associate d to the CMB. In the next section, we pres ent a method, 
based on Hernand ez-Monteagudo & Sunvaev! d2005). that uses the 
caracteristic spectral signature of the ISW-tSZ signal, to separate it 
from the CMB. 
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DETECTION LEVEL OF THE ISW-SZ CROSS 
CORRELATION - IDEAL CASE 



Hernandez-Monteagudo & Sunyaev (2005) have shown that using 
combinations of multifrequency observations of the microwave sky 
can allow to avoid the limit due to the cosmic variance and unveil 
a weak signal whose frequency signature differs from the CMB 
blackbody law. Following this idea, we propose a combination of 
different channels in order to unveil the ISW-tSZ signal, in spite 
of it is dominated by the primary CMB in each channel. Since the 
tSZ vanishes at 217 GHz (in the non-relativistic assumption), the 
channel combination 217 * v, - 217 * 217 (y, # 217 GHz) should 
give an unbiased estimate of the ISW-tSZ,,. cross power spectrum 
in the limit in which we neglect the contaminants such as point 
sources or galactic dust. 

In a first analysis we calculate the signal to noise ratio (S/N) 
of the ISW-tSZ in the ideal case in which the noise is only consti- 
tuted of primary CMB and tSZ autocorrelations. The variance of 

, _ISW-tSZ„. . 

the C t ' is : 

A 2 (C? W " SZ "i) 



,:, 7 ) 2 >-<4 



■ a 217 a 217 ) 2 



which reduces, in the ideal case, to 



a 2 (c; 



i)= [ c ISW-,SZj- + ( c ISW + C CMB + 



N f )(cf 



(23) 



^tSZlh\ 



+ C { ' ). 

(24) 

The kinetic SZ effect has the same spectral dependence as the pri- 
mary CMB therefore it will be subtracted by the channel combi- 
nation we propose. Furthermore, the kinetic SZ contribution, at the 
large scale of interest in this study, is orders of magnitude smaller 
than the one coming from the primary CMB, it thus does not affect 
our estimation of the noise term (Ed. 1 241 ). 

For a full sky coverage, the number of independent modes is 
2( + 1 and can be approximated as / s k y (2^ +1) for a partial sky 
coverage / sky . We therefore write the signal to noise ratio for the 
ISW-tSZ cross correlation at multipole I as 



(21 + l)/ sty [c 



ISW-tSZ 



+ Cf Zlh ) 



[ C ISW-,SZ]- + ( C ISW + C CMB + ^ ( c ,! 

The cumulative S/N up to multipole £ for a full sky survey is 
(2f + l)[c^ w - ,sz f 



1/2 



(25) 



N 



if) 



z 



1/2 



[CfS W - lSZ f + ( C BW + C CMB + ( c ,SZ2h + piSZIh) 

(26) 

In case of partial sky coverage or sky cuts of the Galaxy, the low 
multipoles are not independent anymore, one should adopt a con- 
servative approach starting the summation fro m l m j n ^ n/8 s ^ y , 6 s ^ y 
being the survey size in its smallest dimension dCooray|2 002). Note 
that in the expression above the g 2 , factors describing the frequency 
dependence of the tSZ effect cancel out. This means that in the 
ideal case the choice of the frequency v, does not affect the signal 
to noise level of the ISW-tSZ cross-correlation measurement. Nev- 
ertheless, as discussed in section 5, the level of contaminants will 
depend on the frequency choice which will impact on the signal to 
noise ratio. 

In Fig. [3] we represent the signal to noise ratio contribution 
from each multipole when considering a cosmic variance limited 
full sky experiment. The solid and dotted lines show the results 
obtained when using Eq. J 19b and the Limber approximation of Eq. 
<21t respectively. Most of the contribution to the signal to noise 
ratio comes from multipoles ranging from 5 to 30 with a maximum 
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Figure 3. S/N contribution from each multipole for a cosmic variance lim- 
ited full sky experiment with no contaminants. The solid black is obtained 
without Limber approximation while the dot dashed blue line is obtained 
using the Limber approximation as in Eq. 12 It 




1000 



Figure 4. Cumulative signal to noise ratio, same as Fig. [5] 



contribution at C = 13. This multipole range corresponds to the 
scales in which the ISW-tSZ cross term is expected to reach its 
maximum amplitude and more importantly dominates over the tSZ 
autocorrelation term. 

The results for both the no Limber and Limber approximations 
are presented in Fig. [4] We found that the cumulative S/N is of the 
order of 2.2, w hich is about half the signal to noise ratio found by 
ICooravl |2002). The main difference between the computation of 
the angular powe r spectra in his article and our computation is that 
we introduced the K omatsu &Seliakl (f2002) cluster pressure profile 
(Eq.[l6]l when estima ting C' sw ~' sz a nd Cf z ~ tSZ . It is not clear from 
Eqs. (32) and (33) in lCooravl d2002l) if they account for the cluster 
pressure profile and what model they use. 

As pointed out by re c ent CMB experimen ts dLueker et al.l 
l20ld : iKomatsu et at] l2010l: iDunklev et all |2010|) , the measured 
tSZ power spectrum seems to have a lower amplitude than ex- 
pected from analytical calculations or simulations. Apart from 
point sources residual that could fill the tSZ decrement observed 
at frequencies below 217 GHz, this discrepancy could be ex- 
plained if the erg value used to compute the cluster abundance 
in the spectrum modelisation is overestimated. Nevertheless in 
this case there would be tension between a lower <x 8 value de- 
rived from SZ observations and values derived from primary 
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CMB measurements i Komatsu et alfeofol) or X-ray cluster counts 
dVikhlinin etafl [2009). The disagreement between the observed 
and predicted tSZ power spectrum could also be due to an overes- 
timation of the tSZ effect induced by each individual cluster. The- 
oretical models describing the electronic density and temperature 
profiles (such as KS02), output from hydrodynamical n umerical 
simulations (such as those obtained by Sehg al et"al] d2010h ) as well 
as pre ssure profiles derive d from X-ray observations (such as the 
one in lArnaud et al.l d2009h ) all seem to predict a higher tSZ power 
spectrum than what is effectively observed. Physical processes such 
as supernovae or AGN energy feedback to the intra-cluster medium 
or a non thermal pressure co ntribution could ex plain the lower level 
of observed SZ signal (e.g. Ishaw et alj|2010h . Here we used the 
theoretical KS02 profile to calculate the power spectra (Eq. 1 16b - 
Comparing the effect of various cluster electronic profiles on the 
amplitude or shape of the ISW-tSZ paper is beyond the scope of 
this paper and will be investigated in a subsequent work. 

The main contribution to the cJ SWtSZ variance in Eq. l !24t is the 
[Cf w + Cf MB ) (Cf z2h + Cf zlh ) term, and especially the tSZ con- 
tributions that multiply the CMB noise term. We then propose to 
mask some SZ clusters in order to significantly reduce the noise 
term without decreasing too much the ISW-tSZ signal. 




0.010 
0.001 



Figure 6. The thick lines represent the power spectra after masking clusters 
at z < 0.3 (The color coding is as follows : solid red : ISW-tSZ, dot-dashed 
blue : ISW-ISW, short dashed green : tSZ-tSZ (lhalo term), long dashed 
green : tSZ-tSZ (2halo term)). For comparison the thin lines represent the 
contribution due to all clusters. The dotted line is the primary CMB angular 
power spectrum. 



4 OPTIMISING THE ISW-TSZ DETECTION 

In order to optimise the detection of the ISW-tSZ signal, we need 
to determine which clusters contribute more to the ISW-tSZ signal 
and which clusters are responsible for the tSZ-tSZ shot noise on 
large scales. We therefore compute the halo contributions to the 
CfS, as a function of their redshift, for different £ values. We present 
the results in Fig. [5] On the large scales (several degrees) we are 
interested in, the tSZ signal mainly comes from clusters at a redshift 
lower than 0.3 for the 2-halo term (dot-dashed red line) and 0.03 
for the one halo term (dashed blue line). On those same angular 
scales, the ISW signal (dotted green line) mainly arises from time 
varying potential wells at redshift [0.2; 1]. The cross ISW-tSZ signal 
mainly comes from clusters at a redshift z e [0.06; 0.8] (solid black 
line). We therefore expect that masking the low redshift clusters 
will enhance the signal to noise ratio of the ISW-tSZ. 



4.1 Sharp redshift cut 

We first apply a sharp cut in redshift in order to simulate the effect 
of masking low z clusters. The power spectra for a z = 0.3 thresh- 
old are presented in Fig. [6] As expected, the lhalo term of the tSZ 
contribution is strongly decreased, as is the 2-halo term to a lower 
extent (solid green lines). For multipoles lower than 100, the dom- 
inant tSZ contribution is now the 2-halo term. Since the ISW-tSZ 
signal comes from higher redshift than the tSZ signal, masking has 
a lower impact on its power spectra. The ISW-ISW signal is only 
reduced at very large scales (( < 10) and this relative variation does 
not exceed 25%. 

As shown in Fig. [7] the drastic reduction of the tSZ 1-halo 
term strongly boosts the ISW-tSZ correlation coefficient, which is 
defined as 

wswtsz 

^isw.sz = ^' (27) 

/7^Sw7nSZlh+2h 

This can be easily understood since the tSZ 1-halo, as opposed to 



■g 0.8 - 
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Figure 7. ISWtSZ correlation coefficient. The solid black line is the co- 
efficient after masking z < 0.3 clusters, which drastically reduces the tSZ 
1 halo term at multipoles lower than 100. The dot-dashed blue line is the 
correlation coefficient without any mask. 



the 2-halo term, term does not trace the large scale correlations and 
behaves as noise for the ISW-tSZ detection. 

We then determine what would be the signal to noise ratio 
that could be obtained for different redshift cut thresholds. Fig. [8] 
shows the maximum signal to noise ratio for the ISW-tSZ cross- 
correlation is reached when masking z < 0.3 clusters. A simple cut 
to mask low z clusters allows to increase the S/N from 2.2 to 5.1. 
Should one also mask M > 5 x 10 14 M Q clusters only, the signif- 
icance level would ri se to ~ 6cr. The se detection levels compare 
well to the results of Afshordi (2004[), who claims that it should 
be expected a detection at the level of ~ 7.5<x for an ideal ISW- 
galaxy correlation. According to this work, this S/N is expected to 
be of the order of 5 for an all-sky survey with 10 millions galax- 
ies within < z < 1, provided the redshift systematic errors are 
lower than 0.05 and the systematic an isotropies of the survey do 
not exceed 0.1%. lDouspis et all f2008) showed that an Euclid like 
mission should ideally provide a detection of the ISW at a signifi- 
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Figure 5. Normalized contribution to the CfS as a function of the halos redshift, for different multipoles. The color coding is as follows : dotted green : 
ISW-ISW, solid black : ISW-tSZ, dot-dashed red : tSZ-tSZ (2halos term), dashed blue : tSZ-tSZ (lhalo term). 
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Figure 8. Cumulative signal to noise ratio of the ISW-tSZ detection in 
the ideal case (f s ty = 1, no foregrounds). The dot-dashed curves repre- 
sent the S/N obtained for different cut in redshift : 0.1 (purple), 0.2 (blue), 
0.3 (light blue), 0.4 (green) and 0.5 (red). As a comparison the solid black 
curve reminds the S/N obtained without masking any cluster. The high- 
est solid black curve is the cumulative S/N obtained for a z < 0.3 and 
M > 5 X 10 14 M o cut. 



cance level of 5<r. Note that we obtain similar signal to noise ratio 
with our approach but without the need to use non CMB data. 



4.2 A more realistic approach : using a SZ selection function 

The previously introduced sharp cuts in redshift assumed that red- 
shifts and masses were available for a relevant subset of clusters. 
We now relax this assumption and build a theoretical selection 
function in order to determine which clusters can be detected -and 
thus masked- in a Planck-like survey. 

A galaxy cluster is assumed to be detected through its SZ sig- 
nal if its beam convolved Compton parameter exceeds the confu- 
sion noise and its integrated signal Y is higher than A times the 
instrumental sensitivity simultaneously in the 100, 1 43 and 353 
GHz channels, with A the detection significance level, (Bartelmann 
l200ll;lTaburet et alj|200fl) The selection functions for 1, 3 and 5 a 
detections are represented in Fig. [9] 

As shown in Fig.[T0] masking the clusters detected with a high 
detection significance (5 cr, red line) increases the cumulative sig- 
nal to noise ratio to 3.7. When masking clusters that are detected 
at a 3cr level, the cumulative S/N is enhanced to 4. Reaching this 
level of detection requires to mask faint clusters and thus suppose 
a good understanding of the cluster selection function. 



5 IMPACT OF POINT SOURCES ANG GALACTIC 
RESIDUALS 

In the previous section we considered an idealized scenario in 
which the only signals were the primary CMB, the ISW, the tSZ and 
the instrumental noise. Nevertheless it is of common knowledge 
that radio and infrared galaxies as well as our Galaxy are important 
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Figure 9. Selection functions for 1 (black solid line), 3 (dot-dashed blue 
line) and 5 (triple dot dashed red line) <x detection of SZ clusters for a 
PlanckAike survey. 
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Figure 10. Cumulative signal to noise ratio of the ISW-tSZ detection after 
masking the 5 (dot-dashed red line) and 3 (dot-dashed green line) cr detected 
SZ clusters, black line reminds the result without masking any cluster. 

contributors to the observed signal at submillimeter wavelengths. 
For that, we use state-of-the-art models describing the emission of 
the Milky Way a millimeter wavelenghts, together with the contri- 
bution from radio and submillimeter sources. The frequencies for 
which we built our foreground model are 100, 143, 217 and 353 
GHz, which correspond to the HFI instrument of the Planck mis- 
siorQ. We considered five different contaminants, namely free-free, 
synchrotron and dust emission (coming from our galaxy), and ra- 
dio and infrared extragalactic sources. For the free-free and syn- 
chrotron, we scaled the maps produced by the WMAP team at the 

V and W bands of this experiment (corresponding to 74 and 94 
GHz, respectively) to higher frequencies. We used the maps avail- 
able at the LAMBDA repositorjl 3 ], and computed the effective spec- 
tral index in thermodynamic temperature for each pixel between the 

V and the W bands. These spectral indexes were used to extrapo- 
late the thermodynamic temperature to the higher frequencies un- 
der consideration. The dust emission induced by the Milky Way 
was predicted in our frequency range my means of model 8 of 

2 URL site of HFI: http : //www . planck . fr /heading 1 . html 

3 LAMBDA URL site: http://lambda.gsfc. nasa jjov] 



Finkb einer et alj i 19991 hereafter denoted by SFDfl This model 
predicts the galactic dust emission down to angular scales of ~ 5 
arcmins, which suffices in the context of ISW studies . Final ly, we 
used the point source maps produced bv lSehgal et al .1 12010), who 
modeled the contribution of the extragalactic source population, 
both in the radio and the submillimeter. In that work, a halo cata- 
log resulting from a large cosmological hydrodynamical simulation 
was populated with radio and dusty galaxies, in such a way that 
the spatial distribution of those sources follows the clustering of 
dark matter halos. Both galaxy populations were designed to meet 
different constrai nts obtained from ra dio, infrared and millimeter 
observations, (see lSehgal et al.1 {2010) for details). They were pro- 
jected on high resolution sky maps at different observing frequen- 
cies, which differ from HFI channels' central frequencies for a few 
percent in most of the cases. 

Our proposed approach to unveil the ISW - tSZ correlation 
compares CMB observations obtained at frequencies for which the 
tSZ is non zero with CMB observations observed at 217 GHz. This 
frequency lies in the Wien region of the CMB black body spectrum, 
in which the CMB brightness drops very rapidly for increasing fre- 
quencies and dust contribution (both galactic and extragalactic) be- 
comes dominant. Therefore in a relatively narrow frequency range 
the CMB is surpassed by dust emission, and it is in this frequency 
range where an extrapolation of dust properties (observed at high 
frequencies) down to lower frequencies (where CMB is dominant) 
must be carried out. On the large angular scales of relevance for 
the ISW it is dust in the Milky Way the main source of contamina- 
tion, and its accurate subtraction is actual critical for our purposes. 
An experiment like HFI counts with frequency channels centered 
at 353, 545 and 857 GHz, which probe the regime where dust emis- 
sion is well above the CMB contribution. We shall use those chan- 
nels to correct for dust (both galactic and extragalactic) at lower 
frequencies. Our approach attempts by no means to be exhaustive 
nor systematic, but simply tries to display the degree of accuracy 
required at subtracting dust emission in order to unveil the tSZ - 
ISW cross correlation in a Planck-like experiment . 

We first built a mask that covered those regions where 
the Milky Way emission, both in radio and sub-millimeter was 
stronger. We sorted in intensity (from bigger to smaller values) the 
templates of free-free and synchrotron in the V band (as produced 
by the WMAP team) and the SFD dust template at 353 GHz. Mask- 
ing a given level of emission (for instance, the 25% of brightest 
pixels) in each template yielded two masks that were very similar 
(particularly in the galactic plane), with differences corresponding 
mostly to high latitude clouds being bright either in the radio or 
submillimeter (but not on both). The final mask was the product of 
the two masks built upon the radio and dust templates. The fraction 
of un-covered sky, f^, was then set as a free parameter in the mask 
construction. 

According to the SFD dust templates, one has to take into ac- 
count the spatial variation of the effective spectral index if one is 
to accurately correct for dust emission in the 100 - 217 GHz fre- 
quency range. In these templates, an effective spectral coefficient 
in thermodynamic temperature (defined as the ratio of thermody- 
namic temperatures between two different channels, i.e., 12353, j = 
dTj/ST 35 2GHz) i s correlated with the thermodynamic temperature 
at 353 GHz, as the left panel of Fig.[TT]shows. The curvature at low 



4 The data and the software to build this galactic dust map were down- 
loaded from 

http : //www . astro . princeton . edu/~schlegel/dust/cmb/cmb . html 
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temperatures is a consequence of the grey body law describing the 
dust emission in IR galaxies, and we make use of it when subtract- 
ing the dust emission at low frequencies. In this low temperature 
regime, the effective spectral coefficients from IR galaxies differ 
from the that of the Milky Way, and for this reason a more accurate 
scaling could be obtained by treating the local cirrus component 
separately from the extra-galactic IR one. This can be achieved, in 
high lattitude regions, using HI data that traces the galactic cirrus 
emission, in order to remove their contribution. We sorted pixels 
of HEALPi^ resolution parameter N sjde = 64 outside the mask 
according to their intensity in the dust template at 353 GHz, and 
binned them in groups of length n groups , in each of which a different 
estimate of 0-353,; is estimated in the low frequency channel^ At 
these frequencies the observed signal in pixel n is modeled as 



Galaxy at 317 GHz 



N\h) + 0-353, j(n)M 353 (n) + R J (n), 



(28) 



where M 353 (n) is the dust template at 353 GHz (including both 
galactic and extragalactic emission), R j (n) is the extrapolation from 
the radio template of WMAP's V band, and A"(n) is what we re- 
gard as the noise component at frequency j. Most component sepa- 
ration algorithms attempt to fit for all relevant components (includ- 
ing CMB, radio and dust contributions) at each frequency. Now 
we concern only about the impact of dust residuals, which at these 
frequencies are dominant. We find however that by using a mask 
acting on the 10% brightest pixels in the synchrotron plus free-free 
template, the effect of radio remains always below that due to dust. 
The effective noise N j contains the residuals of a first order CMB 
subtraction, and will be assumed to show no spatial correlations 
(i.e., to be white noise). For a high level of effective noise it will 
be necessary to bin the signal in larger groups (bigger n group!: ) in 
order to measure an average 0-353,; throughout the bin. However, 
this average estimate of the effective spectral coefficient in the bin 
will produce an estimate of the dust contribution at frequency j and 
pixel n that will not correspond exactly to the exact value, and this 
mismatch will be more relevant the wider the bins are (i.e., the big- 
ger ngroups ls )- On the other hand, if noise were negligible one could 



make n„ 



1 and measure the dust contribution very accurately 



at each sky position n. 

The residuals from the dust subtraction can be then written as 

ST^Jn) = THn) - a 353 , jM 353 (h), (29) 

where the 0-353 ; estimate is computed after minimizing 

X 2 = X" - & 353, jM 353 )^ (V - 0353. jM 353 ) h (30) 
M 

and the indexes i and Z are running from 1, n gr0 ups in the bin to which 
the pixel n belongs. The matrix N -1 corresponds to the noise in- 
verse covariance matrix, which is diagonal if A' is poissonian noise. 
The estimate and error associated to 0-353 j are given by 



tt353. 



2y M 353 NjJ M 353 ' 

1 



(3D 
(32) 



Z u M 353 NjjM 353 ' 

The residuals given by Eq(29] will bias the tSZ - ISW corre- 
lation estimates, and also increase their errors. The former is given 



5 HEALPix's URL site: http : //www . healpix . org 

6 Unless otherwise specified, this is the pixelization resolution used in our 
analyses 





Residuals at 217 GHz 
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Figure 12. Top panel: Total galactic emission outside the mask covering 
40% of the sky at 217 GHz in thermodynamic temperature units. Bottom 
panel: Residuals outside the mask at 217 GHz after using a perfect dust 
tracer at 353 GHz and assuming a white noise level of 80 pK at scales of 5 
arcmins. 



by the second term in right hand side of this equation: 
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(33) 

where the multipoles of the residuals at channel j (Eq.|29t are given 
by the a!'™ s -s. The impact of residuals on the dispersion of the tSZ 
- ISW correlation estimates is modeled as in Eql24l 

= (cf w -' sz 'f + c« s (cf z " "' + cf z * 2 " + cf 217 - res ), 

(34) 

the power spectrum multipole of the residual differ- 
ence map ST 3 res - 8T™ and where we have assumed that the power 
spectrum of residuals at 217 GHz is much smaller than the CMB 



A[ C 'r- ,sz '] 



with C 



7-217, 

; 



power spectrum (Cj 



«c Cf MB ). As before, the S/N for a given 

- tSZ cross correlation 
bin size n gmup and Gaussian 



multipole will be given by (S/N)"f = (C ; 

We computed the total S/N for the ISW 
for different levels of sky coverage f sky 
white noise amplitude cr N = ^(Na). We considered the case where 
all clusters below z < 0.3 are masked out, since it is the one giv- 
ing rise to highest S/N in the ideal (contamination free) case, and 
a low frequency of 143 GHz. We first assumed a noise amplitude 
of o"jv = 80 pK in square pixels of 5 arcmin size, and checked for 
the resulting S/N after binning in groups of varying lenght n groups . 
This noise amplitude is above the predicted values for HFI chan- 
nels 100, 143 and 217 GHz, for which values of 14.2, 9. and 14.6 
p: K in pixels of size 5 arcmins are expected. One would guess, 
however, that any component separation algorithm (required to re- 
move an estimation of the CMB in each channel) would sensibly 
increase the overall noise level. We found that after binning pixels 
in groups of size n groups = 1000 and masking 40% of the sky, resid- 
uals would drop below the ISW - tSZ cross correlation amplitude 
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Figure 11. Left panel: Effective spectral coefficient between 353 GHz, on the one hand, and 143 GHz (green points) and 217 GHz (blue points). Middle panel: 
The expected ISW - tSZ cross power spectrum amplitude (at 143 GHz, after masking all clusters below z < 0.3) given by the solid black line is compared to 
the residuals present in the 143 GHz (green circles) and 217 GHz (blue circles) channels, respectively. The effective bias in the ISW - tSZ cross spectra are 
displayed by red circles, while the black ones display the increase in the error budget (last term in Eq |341 . Units are expressed in Di = /(/ + l)C//(2;r). Right 
panel: Cumulative S/N of the ISW - galaxy cross power spectrum below a given multipole /. 



for the 143 and 217 GHz channels. In the middle panel of Figll II 
green circles display the amplitude of dust residuals at 143 GHz, 
while those at 217 GHz are given by blue circles. The effective bias 
on the ISW - tSZ cross correlation (which is given by the solid 
thick line) are shown by red circles. The extra (last) term in the er- 
ror computed in Eq[34]is displayed by black circles. In all cases we 
are showing pseudo-powei spectrum multipoles. In this configura- 
tion, the total S/N achieved was approximately 2.6 ( as shown by 
the rig ht panel of Figll lb. When comparing to Fig. 6 of lLeach et al.l 
(2008), we see that our level of residuals is comparable to those 
obtained after using the whole set of Planck frequencies for com- 
ponent separation. The accuracy of our simple cleaning procedure 
is displayed by Fig|T2j the top panel shows the total contaminant 
emission at 217 GHz, while the bottom one displays the residuals 
after the cleaning procedure described above was applied. 

We also explored how S/N depends on cr N and n grmps , see 
Fig. [T3] The top panel shows how S/N varies with cr N for differ- 
ent values of n grmlps : black, red, green and blue lines correspond 
to n groups =\,\Q, 100 and 1000, respectively. In the bottom panel, 
those same colors correspond to cr N = 1, 10, 50 and 100 /iK re- 
spectively. The top panel displays how S/N is degraded as the noise 
level is increased, being this degradation more important for low 
values of n gmups for the same value of a N . However, in the cases 
of low noise, the reconstruction of dust is more accurate for small 
Kgroups an d higher values of S/N are reached in such cases. The bot- 
tom panel, instead, shows how increasing the bin width (n grm<ps ) 
improves the S/N only in the regime where noise is dominating: 
for low noise levels binning only degrades the quality of the sub- 
traction. Note that in ideal conditions (n groups = 1, cr N = //K), 
S/N^ 3.9 for f sky =0.6 considered here. 



6 DISCUSSION AND CONCLUSIONS 

The problem of component separation in multifrequency CMB 
observations for experime nts like Planck has been a subject of 
active investigation (e.g., IStomporetaIll2009l ; iLeach et alj|2008l ; 



Eriks en et alj2008l ; lAumont & Maci'as-Perezl2007tlStolvarov et all 
2005). The different contaminants (either galactic and extragalac- 
tic) show generally a different spectral dependence when compared 
to the CMB. On the low frequency side, foreground powerful in 
radio wavelenghts fall steeply and become th e sub-dominant fore - 
grounds at frequencies above v ~ 70 GHz, jBennett et ai]|2003h . 
Above this frequency, the presence of dust absorbing UV radia- 
tion and re-emitting it in the sub-millimeter and millimeter range 
constitutes the dominant contaminant. An experiment like Planck, 
with four high-angular resolution channels covering the 217 - 857 
GHz frequency range, should provide an accurate description of 
this foreground. With these data in hand, currently existing models 
describing the physics of dust emission should be improved further 
and accurate extrapolations to lower frequencies should be enabled. 

In our toy model describing the impact of the contaminants 
we assumed that the 353 GHz channel was a perfect (CMB - free) 
dust tracer, and that the signal at low frequencies (100 or 143 GHz) 
was either due to dust or a white Gaussian noise signal. These as- 
sumptions may be overly optimistic (in regard to the 100 & 217 
GHz channels themselves), but however a generic combination of 
all nine channels (ranging from 30 GHz up to 857 GHz) are ex- 
pected to provide an estimation of each of the components (CMB + 
foregrounds) wh ose residuals are e xpected to lie close to the level 
of ~ 0.1 (pK) 2 jLeach et al.l [20081) . This is already the accuracy 
ballpark that our simple analysis proved to define the regime of 
detectability of the tSZ - ISW cross correlation, and there may 
still be room for a more optimized channel combination oriented 
to unveil the particular tSZ-ISW cross correlation. The use of high 
galactic latitude HI maps as tracers of galactic cirrus could allow 
to lower the impact of those residuals at t he level of a percent in 
Ci, (G . Lagache, private communication). iFernandez-Conde et al.l 
(2008) computed the cirrus power spectrum at different frequen- 
cies and different HI column densities. They found, in the case of 
fields that have a very low level of dust contamination, that the 
cirrus power spectrum at 217 GHz is of the order of 5 (juKrj) 2 at 
£ = 10 in units of /(/ + l)C//(27r). A one percent residual of the cir- 
rus emission is thus smaller than 0.05 (/jKrj) 2 , i.e., 0.4 (//Kcmb) 2 - 
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Figure 13. Top panel: S/N versus residual effective noise for different val- 
ues of ngroups'' black, red, green and blue lines correspond to n groups =\,\0, 
100 and 1000, respectively. Bottom panel: S/N versus bin width n gmups : 
black, red, green and blue lines correspond to cr/v = 1, 10, 50 and 100 /iK 
respectively. 



Current foreground residual esti mates based upon t he works of 
iFernandez-Conde et all J2008I) and lLeach et al.l ([2008) suggest that 
at the frequencies of interest (100 - 217 GHz) the contaminant 
residuals remain a factor of a few above our requirements. How 
much room there is for improvement below those limits is some- 
thing yet to be estimated from real data. 

Planck data at high frequencies provide a more profound 
knowledge of the dust properties in both our galaxy and extragalac- 
tic sources, together with the mechanisms involving its emission 
in the sub-millimeter range. It is nevertheless important to bear 
in mind that, since the different frequencies sample the infrared 
galaxy populations at diffe rent redshifts and the galaxy linear bia s 
evolves with redshift (e.g. lLagache et ai]|2007l ; Iviero et alj|2009h . 



using high frequency maps so as to clean the 143 and 217 GHz 
could potentially degrade the residual level we obtained in Section 
5 with our template fitting method. This issue is under current in- 
vestigation within the Planck collaboration. 

Even in the worst scenario in which foreground residuals are 
too high and complicated to prevent the detection of the tSZ - ISW 
cross correlation, the upper limits to be imposed on it are of cosmo- 
logical relevance, since it would constrain cosmological parameters 
like cr 8 , Q. m or S1 A . 

We have shown that the tSZ - ISW cross correlation con- 
stitutes a CMB contained test for Dark Energy. The peculiar fre- 
quency dependence of the tSZ effect and the availability of multi- 
frequency all sky CMB observations provided by the experiment 
Planck should enable an estimation of this cross correlation pro- 
vided the hot gas is a fair tracer of the potential wells during the 
cosmological epochs where the ISW is active. Our theoretical study 
shows that the Poisson/shot noise introduced by the modest num- 
ber of very massive, very bright in tSZ galaxy clusters can be at- 
tenuated by masking out those tSZ sources below redshift z < 0.3. 
In the absence of a massive galaxy cluster catalog below that red- 
shift, it would suffice to excise from the analysis those tSZ clusters 
clearly detected by Planck in order to achieve S/N of the order of 
3.9 (f s ky = 1). This tSZ - ISW cross correlation detection would 
not require the use of deep-in-redshift and wide-in-angle galaxy 
surveys, but only the combination of different frequency CMB ob- 
servations. This would hence provide a different approach for ISW 
detection with different systematics to other attempts based upon 
CMB - galaxy survey cross correlations. If foreground residuals 
are kept at or below the ~ 0.04 QiK) 2 level (in /(/ + l)C//(27r) units 
at / ~ 10) in the frequency range 100-217 GHz, then the tSZ - 
ISW correlation should provide a valid and independent test for the 
impact of Dark Energy on the growth of structure and the evolution 
of large angle CMB temperature anisotropics. 
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